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ABSTRACT 



We suggest a simple analytic approximation for magnitudes and hence distances 
of Blue Horizontal Branch (BHB) stars in Sloan colours. Precedent formulations do 
not offer a simple closed formula, nor do they cover the full dependences, e.g. on metal- 
licity. Further, using BHB star samples from the Sloan Digital Sky Survey (SDSS) we 
validate our distance calibration directly on field stars instead of globular clusters and 
assess the performance of other available distance calibrations. Our method to statis- 
tically measure distances is sufficiently accurate to measure the colour and metallicity 
dependence on our sample and can be applied to other sets of stars or filters. 



Key words: 

Galaxy: halo - 



stars: horizontal branch - stars: distances 



1 INTRODUCTION 

Blue Horizontal Branch (BHB) stars a re among the most 
popular tracers o f the stellar halo (e.g. iPreston et al . 1991:; 
Isirko et al.l |2004 ). They are bright enough to cover dis- 
tances larger than 5 kpc and span a relatively narrow range 
in a bsolute visible magni t ude, if one avoids the blue hook 
(see 



Whitnev et al] 1 19981 : iD'Cruz et all I2OO0I : iBrown et all 



I2OOI . for a definition), where the luminosities of the blue- 
most BHB stars are greatly diminished. Further they have 
very blue colours, reducing the risk of being confused with 
other types of stars and are globally available - in contrast 
to RR Lyrae stars, which, in addition to spectroscopy, re- 
quire a good cov erage in time to resolve their oscillations 
l|Sesar et al.ll2010l ). 

By far the largest sample of BHB stars comes from 
the Sloan Digital Sky Survey (SDSS. lEisenstein et al]|2011^ . 
Some previous studies attempted to provide a magnitude 
calibration for BHB stars in the ugriz filter system of 
SDSS, based on a combination of theore tical expectations 
and calibrating on globular clusters (e.g. ISirko et al.ll200j : 
IXue et al]|2008l ). Yet those calibrations have not been pro- 
vided in a user-friendly form, such that many current 
studies (e.g. iDeason et al.l l2011al : iDe Propris et all I2OI0I : 
iNiederste-Ostholt et al.ll201(]| )" still use the approximation of 
a constant ab solute magnitude Mg =0.7, originally derived 
in lYannv et'al . (2000) by comparing RR Lyrae and BHB 
stars in the globular cluster Pal 5. 

More impo rtantly, the classic calibrations (like 
ISirko et al.|[2004 ) rely on globular clusters, bearing the usual 
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uncertainties in their distance moduli and reddenings. Fur- 
ther we know that (a) the Horizontal Branch (HB) stars 
of clusters vary between clusters for reasons that are in- 
completely understood (e.g. clusters of the same metallicity 
can either be very red or blue: specifically metal-poor clus- 
ters with high central density present bluer HB morphologies, 
fSuntzeff et al..,199L : Buonanno 1993. ) and (b) there are sys- 
tematic differences between globular-cluster and field halo 
stars (e.g. the binary fraction among extreme HB stars is an 
order of magnitude l ower in globular clusters than among 
field stars. llfan|[2Cioi ). 

Given these problems, it is crucial to create a formula 
that fulfils the following criteria: 

• it covers the principal dependences of the HB magni- 
tudes (i.e. colour and metallicity) 

• it is provided as a simple analytical formula that gives 
the absolute magnitude from observables; 

• it is validated by /calibrated on field stars. 

The aim of this paper is to construct such a calibration 
in the ugriz filter system of SDSS and assess its performance 
w.r.t. currently available calibrations by testing it on BHB 
field stars. At first glance the last point seems out of reach, 
considering that there are virtually no reliable Hipparcos 
parallaxes for BHB stars and especially none for stars within 
the magnitude range of the SDSS. However, the sample size 
of the SDSS enables us to obtain statistical distance esti- 
mates for fi eld BHB stars v i a the V-W correlation method 
laid out bv lSchonrich et all (|2012l ). Even though this exer- 
cise is performed on SDSS/SEGUE, our method is general 
and can be applied to any stellar sample of sufficient size 
and quality in any filter set. 
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Models of stellar evolution predict the locus of the HB. 
Hence, our empirical calibration will also provide a test 
of these models. We find that our cali bration agrees quite 
close ly with BASTI isochrones (jPietrinferni et all |2004 
|2006| ). and makes stars signifi cantly more luminous t han 
Dartmouth isochrones predict (jPotter et al ] |2007l . l2008h . 



2 GENERAL METHOD 

We will first identify a general method to assess the ac- 
curacy of any distance calibration and then formulate one 
that fulfils the requirements outlined in the introduction. In 
principle one could use the distance measurement in the fol- 
lowing subsection directly to produce a distance calibration 
that solely relies on the measured data. As we will see on 
the BHB star d ataset our distance calib r ation derived from 
BASTI models l|Pietrinferni et al.l |2004 120061 ') is a signifi- 
cant improvement w.r.t. to precedent approaches. Since it 
is within the measurement errors indistinguishable from the 
optimum fit to the field star data, we stick with the formula 
directly derived from theoretical models. 

2.1 How to assess a distance calibration 

ISchonrich eraH (|2012l ). hereafter SBA12, investigated how 
distance errors introduce correlations between velocity com- 
ponents. They consider a distance fractional error / such 
that the measured average distance is s' — {1 + f)s, where s 
is the true distance. The observed {U, V, W) velocity compo- 
nents can then be expressed as functions of both the real ve- 
locities and of /, which introduces correlations between the 
velocity components. The main method of SBA12 is not ap- 
plicable for BHB stars, since the (uncertain) proper-motion 
errors enter the correlation multiplied with the square of 
the distance. In practice, with current proper motions the 
method is applicable to samples of stars with s < 5kpc. 

However, we can directly use the linear estimator given 
in their equations 42-45. This estimator uses the fact that 
halo stars have a large mean heliocentric azimuthal (V) ve- 
locity because the halo does not have any significant rotation 
so we see the full reflex motion of the Sun. Any mean dis- 
tance error then translates this mean motion into a mean 
motion in the radial U and vertical W components of he- 
liocentric velocity depending on galactic angles according to 
the matrix elements provided by SBA12. For the correlation 
between mean rotational motion and vertical heliocentric 
mean velocity we have (equation 45 of SBA12): 

{W) + Wq=c- Twv + e (1) 

where c = f{V)/{l + f{Tvv}) is the slope of the linear best 
fit between {W) + Wq, while Tyv = 1 - cos^bsin^/ and 
Twv = — ^sin2fesinZ are the angular terms through which 
distance errors introduce correlations between the apparent 
components of velocity (see Table 1 of SBA12). In particular 
a correlation between W and Twv implies that the stellar 
halo is crowding towards the disk on one side of the Galaxy 
(w.r.t. the Sun-GC line) and dispersing away from the disk 
on the other side (the "falling sky " effect, in the nomencla- 
ture of SBA12). 

In practice, we perform a linear least squares fit of 
(W) + Wq versus Twv to estimate the correlation term 



c. Then, we iteratively multiply all distances by correc- 
tion factors until this correlation on the sky disappears. 
The fractional average distance error is then just given by 
the correction factor that makes this correlation zero and 
the uncertainty on this estimate directly translates into the 
error bars of c. As the method is linear, random proper- 
motion errors are irrelevant, apart from increasing the error 
bars by expanding the measured dispersion in W. We are, 
however vulnerable to systematics in the proper motions. 
Hence we use th e prop er-motion correc t ion w orked out by 
ISchonrichI (|2012l ) on the lSchneider et al.1 l|2010l) quasar sam- 
ple in Galactic angles {l,b). 

2.1.1 Internal accuracy of the method 

We assess the internal accuracy of the method by testing 
its ability to recover the mean fractional distance error on 
a range of control cases. In each of these, we introduce an 
input bias in the distances and proper motions of the stars 
in the XI 1 sample so that the average distance offset is /in: 

(*^)^/.„ , (2) 

where di = diCnra,n{Qg{li,bi, di) is the perturbed distance 
of the i-th star, g{li,bi, di) accounts for the angular position 
of the star in the sky, ran((^) is a random number generator, 
and c„ makes sure that ((2} is satisfied. 

We select eight different values of /in in the range 
[-0.2,0.2] (i.e. ±20%). Using the V-W estimator to es- 
timate the mean fractional distance error on the sample, 
we recover the true bias within one standard deviation: 
/in - fv-w = 0.0162 ± 0.0187. 

2.2 Calibration construction 

Stellar evolution models predict a steep decline of the in- 
trinsic magnitude of BHB stars towards bluer colours (e.g. 
BASTI, Dartmouth isochrones) in concordance with pho- 
tometry of globular clusters. Further, they predict metal- 
poor BHB stars to be systematically brighter than their 
metal-rich counterparts. The ideal calibration would encom- 
pass the main dependences with the simplest possible form: 
in particular, the metallicity dependence, modelled best by 
a polynomial for the increasing trend in metallicity, and the 
steep decline on the blue side, which we model by an expo- 
nential in {g — r)o colour: 

M,{{g - r)o, [Fe/H]) = aoexp(ai(3 - r),,) + p([Fe/H]) (3) 

where p([Fe/H]) G R„[[Fe/H]]. Parameters of ([Sjl are empir- 
ically adjusted to mimic both the trends seen in metal-poor 
clusters and the ones predicted by the isochrones: among 
the ones available the BASTI isochrones drive our calibra- 
tion towards a better performance on field stars w.r.t. for 
example the Dartmouth ones (see ijS.ip . One must account 
for photometric and reddening uncertainties though, which 
make it likely that a star with assigned blue-hook colours 
will conceivably be still relatively luminous: thus the pro- 
posed calibration shall bend a bit less steeply towards the 
blue hook than the isochrones. 

Our calibration will not be a fit to the data, but rather 
it will be validated by data. The crucial difference is that a 
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priori we do not know if a data sample is perfectly unbiased 
and how eventual systematics can bias the parameters of 
((3)). On the other hand a validation method can rely on the 
limited information content, whi ch proves to be sy stematics- 
free. In the BHB stars sami aleof IXue et Zl (|201li ) for exam- 



ple, iFermani fc SchonrichI }2Q1± flag a peculiar kinematics 
associated with metal-poor high stars: these stars affect 
the correlation between the V component of velocity and 
the angle term sin(Z) sin(6) cos(6), but have negligible effect 
on the relationship between W component of velocity and 
the same angle term arising from the rest of the sample. 



3 SAMPLE SELECTION 

SDSS and its sequel Sloan Extension for Galactic Under- 
standing and Exploration (SEGUE) obtained an unprece- 
dented photometric and spectroscopic mapping of the Milky 
Way galaxy, allowing the study of nearby dwarf popula- 
tions as well as giant star populations up to 100 kpc from 
the Sun (jYannv et al.ll2009l ). The SEGUE Stellar Parameter 
Pipeline provides line index measurements for the spectra 
and the three primary stellar parameters: effective temper- 
ature (TrfF}2_surfece_graviti;_(l^ ([Fe/H]) 
(cfr. iLee et al.ir2008al lbl: lAUende Prieto et alJ lioOSl. Aher- 
native methods have been developed to estimate the above 
parameter s: we will adop t the estimates produced by the 
method of IWilhelm et all (|l999l ) as the latter was particu- 
larly designed for hot stars. 

We build a photometrically selected sample of 
BHB stars dr awn from the 9th Data Release (DR9, 
lAhn fc al.llioT^ '): we query for BIfB candidates and filter in 
colour (dereddened) and stellar parameters in order to min- 
imize the risk of contamination from ot her luminosit y type s 
acc ording to the war nings laid out in lYannv et all (|2000l ') 
and lSirko et al.l (|2004l '). Specifically our selection criteria are 

5< 18 

2 < log(g) < 3.5 
0.8 < {u~g)o < 1.4 
-0.4 < (g-r)o < 
7250 < Teff/K < 9700 
^Sgr 



(4) 



where the last criterion removes 786 stars or respectively 
30% of the sample within a polygon on the sky that encloses 
Sagittarius, the coor dinates of which were kindly provided to 
us by the authors of 'Peaso n et al] (|2011al ) and are reported 
in iFermani fc Scho nrich ( 2 0121 ). These yield a sample of ^ 
2600 stars. A d etailed discussion on relaxin g these criteria 
can be found in IFermani fc SchonrichI (|2012l ). 

We also consid er the spectros copically selected 
sample compiled bv lXue et al.l (|201ll ) and kindly provided 
to us by Xiangxiang Xue on req uest: this consists of ~ 4000 
objects drawn from SDSS DR 8 (lEisenstein et al.ll201ll ) that 
pass both the colour cuts of lYannv et al. ( 200C ) and the 
two B almer-line profile cuts already d escribed in Xue et al] 
l|2008[) , but here slightly relaxed (see IXue et al.ll201ll ). We 
retrieve these objects in DR9 and adopt the latter values for 
both their physical parameters and their astrometry. 

As both samples cover a distance range out to ~ 50 kpc, 
the error in proper motion measurements for single stars is 
approximately as large as the true signal from solar reflex 



motion and velocity dispersion. However, by large number 
statistics the proper motions still bear valuable information. 
While the statistical uncertainty simply sets our formal er- 
ror bars, we have to be concerned with the systematic terms 
that are roughly a factor of 20 smaller than the noise (see 
§4 of lFermani fc Schonrichll2012l ). To cope with the system- 
atics, we will us e the p roper-motion co r rectio n derived by 
ISchonrichI (^201?) on the lSchneider et al.l l|2010l ) quasar sam- 
ple in Galactic angles (Z, b) and note that this correction 
shrinks the statistically predicted mean fractional distance 
error by 0.02 (see §3.1.2). 

Hereafter, we shall refer to the photometric sample se- 
lected via (|4]) as FS12 and to the spectroscopic one as Xll. 
We will assess the distance calibrations on these two sam- 
ples, though we will mainly rely on Xll as it allows for better 
statistics. 



4 DISTANCE CALIBRATION FOR THE BHB 



The end product of the study outlined in i|2.2l is 



M,((5-r)o,[Fe/H]) 



0.0075 exp(-14.0(5 - r)o) -I- (5) 
+0.04 ([Fe/H] + 3.5)^ + 0.25, 

which we show in Fig. [1] for two metallicities: [Fe/H] — 
-2.25 (thick blue line) and [Fe/H] = -1.25 (thick red 
line). Also shown in F i g. [T1 are the BASTI isochrones 
(|Pietrinferni et al.l |2004| . 12009 ). The discrepancy between 
the constant-luminosity approximation (horizontal light 
blue line in Fig. [TJ and both the isochrones and our cali- 
bration exceeds the frequently adopted uncertainty band of 
EA/g ~ ±0.18 (light blue dotted lines in Fig. [T]Q. 

Equation ([5} does not depend on the u — g colour as in 
our tests we could not detect any significant systematics for 
it once the g — r and [Fe/H] dependences had been taken into 
account. The ideal calibration would account for more de- 
tailed chemical dependences like alpha enhancement though 
and particularly the abundance of Helium, which changes 
the absolute magnitude of a star. However, this information 
is generally not available for BHB stars as even though they 
are very hot their Helium ab undance can be very low due 
to depletion mechanisms (e.g. lFaulkner|[l967l ). 

4.1 Validation on field stars 

We validate our calibration on our two samples of halo field 
BHB stars: on both we obtain a satisfactory accuracy: 

/psi2 = 0.01 ± 0.04 and /xii = 0.02 ± 0.03. 

Given that the fractional average distance error from ([5} is 
consistent with zero within the statistical uncertainties, our 
calibration can be regarded as the first direct measurement 
of the magnitude of field BHB stars in SPSS-SE GUE. 

Use of the proper motiorlfl correction of ISchonrichI 

^ An uncertainty of ±0.18 is commonly associated with the cal- 
ibration of the horiz ontal branch, due to neglect of the tem- 
perature dependence l lSirko et al.ll2004l) . or (±0.15, 0.2) associ- 
ated with the constant absolute magnitude approximation (e.g. 
iDeason et~alll2011al : iDe Propris et al.ll2010l) 

^ We checked that the line-of-sight velocity correction given in 
the same paper has virtually no effect on our distance estimates. 
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BASTI, 2 = 0.001 
BASTI, Z = 0.0001 
Mg(g-r,-1.25) 
M,(g-r, -2.25) 
0.7 



-0.2 -0.15 

(g-r)o 

Figure 1. Luminosity as a function of colour for the Blue Hor- 
izontal Branch at different metallicities. Thin lines represent the 
predi ctions from the BASTI isochrones l lPietrinferni et aLll2004l . 
l200d ) as do the dots that are associated with a mass spac- 
ing of 3 ■ 10~®Mq. Thick lines depict our colour and metallic- 
ity dependent approximation for Mg. Blue lines refers to metal- 
poor stars ([Fe/H] G] — 3, —2]) and red lines to metal-rich stars 
([Fe/H] G] — 2, 0]). The light blue line depicts the popular approxi- 
mation of constant Mg =0.7 and light blue dotted lines represent 
an uncertainty of ±0.18. 



reduced the average systematic distance error from 
/ = -0.03 ± 0.04 to / = 0.01 ± 0.04 for the FS12 sample 
and from / = -0.04 ± 0.03 to / = 0.02 ± 0.03 for Xll using 
our calibration (O for the distance estimate. Hence, the nec- 
essary correction affects the average distance calibration by 
about 5% and hardly affects the relative distances between 
single bins. 

In general the statistics are robust against outliers in 
the sense that when we remove stars with extreme proper- 
motion (i.e. tangential velocities above 900kms~^: less than 
3% in FS12 and than 6% in Xll), / does not change for FS12 
and changes by less than 0.01 for Xll. The inclusion of Sgr 
makes very little difference (A/ ~ 0.02) to our estimates. 
Indeed, adding the Sagittarius region to the samples, we find 
/ = 0.04 ± 0.03 on Xll and a similar difference on FS12: 
/ = 0.02 ± 0.04. 



5 DISCUSSION OF OTHER 

MAGNITUDE/DISTANCE ASSIGNMENTS 

Using the same method that allowed us to validate our cal- 
ibration against field halo BHB stars (see tj2.1|l . we can di- 
rectly assess other distance calibrations. 

5.1 BASTI versus Dartmouth isochrones 

In Fig. [2] we plot the discrepancy between the abso- 
lute magnitude predicted by the Dartmouth isochrones 
and our calibration as a function of colour for four 
synthetic populations of ~ 2700 HB stars, each 
with different metallicity. The synthetic populations 
were generated with the Dartmouth web tool (see 
|http ://stellar. dartmouth.edu/models/shb. html [) and input 
parameters ( [a /Fe] , max M,{M), om ) = (0.4,4.0,0.5,0.1), 
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Figure 2. Discrepancy between the absolute magnitude pre- 
dicted by the Dartmouth Isochrones and the one estimated via 
equation [5] as a function of colour for four synthetic populations 
of ~ 2700 HB stars each with different metallicity. Upper panel: 
metal poor stars; bottom panel: metal rich stars. Error bars show 



where M is mass in solar masses. In the colour range g — r£ 
[—0.25, 0] (where our sample stars are concentrated) the off- 
set is systematic and increases with metallicity: AAfg ~ 0.1 
for metal-poor objects and AMg ~ 0.3 — 0.4 for metal-rich 
objectsif] Given that we achieve a fractional average distance 
error consistent with zero well within la on our sample of 
field stars, estimating the stars' absolute magnitudes with 
the Dartmouth isochrones corresponds to underestimating 
distances by up to 4 % for metal-poor stars and by up to 
17% for metal-poor ones (see eg. l5.5.l1 belowl. 

In conclusion, the close agreement between our formula 
and the BASTI isochrones (see Fig. [T} implies that at least 
in the ugriz-system, they are superior to the Dartmouth 
isochrones. 



5.2 Sirko et al. 2004 



ISirko et al.l (|2004h estimated the absolute magnitudes of 
BHB stars as functions of colour for two metallicities 
([M/H] = -1 and [M/H] = -2) by integrating the Ku- 
rucz model atmosphere^ and then finding the point on the 
theoretical track that was closest to the observed star in 
colour space. They validated their calibration with two glob- 
ular clusters: NGC 2419 and Pal 5. With literature distance 
moduli for the clusters they found that the isochrones were 
too bright by 0.13 mag for NGC 2419 and and too faint by 
0.2 mag for Pal 5. They considered that these offsets could 
be accounted for by errors in the adopted distance moduli. 

When we compare our calibration with theirs, it turns 
out that their relation is significa ntly fainter than ours: 
/sirko = -0.07 ±0.04 on Xll. Hence ISirko et al.l (120041 ) pre- 
dict shorter distances than our formula, which on our data 
produces a fractional average distance error consistent with 
zero well w ithin 1 a (se e §3.1) . Further application of the cal- 
ibration of ISirko et al] l|2004h (from their Table 2) requues 

^ Varying \a/Fe\ by ±0.2 does not alter our conclusions. 
* see |http:// www . kurucz . harvard . edu | 
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Figure 3. Average fractional distance error as function of metal- 
licity on Xll sample, using botli the calibration of Xue et al. 
(2011, red line) and equation [5] (blue line). 

several approximations: one needs to bin in effective temper- 
ature, surface gravity, colour and metallicity. Consequently 
the calibration is undefined for more than half of the stars 
in the sample 

Our calibration is an analytic formula and not a com- 
plicated fitting process to the theoretical models. Since it 
produces significantly smaller distance errors than the cali- 
bration of Sirko ct al. (2004 ). it is not only more convenient, 
but also more accurate. 

5.3 Xue et al. 2011 

The catalogue of BHB stars provided by IXue et al.l (|201ll ) 
includes dis t ances estimated with a method very similar to 
ISirko et al.l (|2004| ): i.e. the authors determine the point on 
the theoretical tracks that is closest to the observed star 
in u — g, g — r space and from the former estimate Mg. 
On Xll, i.e. the catalogue of lXue et al] l|201ll ) with the Sgr 
region removed, their distances appear to be weakly biased 
towards distance underestimates: / = —0.05 ± 0.03. 

Marginalizing over metallicity, the fit to the colour de- 
pendence is comparable between their calibration and ours, 
within the formal errors. However, our formula ((Sjl proves 
to be significantly more accurate at low metallicity when 
we marginalize over colour: Fig. [3] shows the different trend 
in metallicity of / for the two distance calibrations on Xll 
sample. In the range [Fe/H] G [—3, —1.8] the calibration of 

^ To avoid binning the data, we comp ared the two form ulas 
at fixed metallicity. We use [Fe/H] while ISirko et al.l ||2004|) use 
[M/H], but this has no relevance in our discussion. When equating 
[Fe/H] and [M/H] and defining AMg = M9(Sirko) - Mg(FS12), 
we find AMg 6 [0.07, 0.35] for metal-rich stars and AMg e 
[0.15, 0.41] for metal-poor stars. When we assume [M / H] = 
[Fe/H]+0.3, which roughly corresponds to the'Sal aris et al] lll993l ) 
approximation [M/H] = [Fe/H] + log(0.638 exp([a/Fe]) +0.362) 
with [a/Fe] ~ 0.4, the estimates ab ove are unchanged . In conclu- 
sion, given their fainter magnitude. ISirko et al.l | |2004| ) calibration 
will produce shorter distances w.r.t. JSjl by a A/ of at least —0.03 
and up to —0.17. Thus, it is reasonable to expect that on our 
samples their calibration would produce a distance error in the 
range / £ [—0.15, —0.02], regardless of how data are binned. 



IXue et alj (|201ll ) is associated with an average fractional 
distance error of / ~ —0.12 ±0.04, while equation ((5| yields 
/ ~ -0.02 ± 0.04. 

5.4 Deason et al. 2011b 

iDeason et alll|2011bl ) propose a calibra tion of BHB sta rs de- 
rived from 10 star clusters published in lAn et alj l|2008l ) cov- 
ering a range in metallicity of [Fe/H] G [—2.3, —1.3]. They 
detect a colour, but no obvious metallicity dependence: given 
the considerable uncertainties in the adopted distance mod- 
uli and reddenings (0.20 ± 0.04 mag or ~ 10 — 15%, e.g. 
iGratton et al.|[T997l ) and the systematic differences between 
clusters and field halo stars though, this result is not enough 
to rule out a detectable metallicity dependence on field BHB 
stars. The calibration of|Deason et al.1 l|2011bl ) is systemati- 
cally fainter for metal-poor stars and brighter for metal-rich 
stars than ours and leads to a worse performance on the 
FS12 sample. 

On the Xll sample we estimate a global fractional aver- 
age distance error of / = — 0.02±0.03 using their calibration 
(/ = 0.02±0.04 for metal-rich st ars and / = -0.08±0 .06 for 
metal-poor stars), while on FS12 lDeason et al] (|2011lJ ) gives 
/ = -0.03 ± 0.04 (/ = 0.01 ± 0.05 for metal-rich stars and 
/ — —0.05 ± 0.06 for metal-poor stars). When investigating 
the detailed morphology of the halo, this metallicity-related 
bias is both significant in magnitude and relevant for the re- 
search goal. E.g. substructures can be blurred and distorted 
by several kpc. In particular, the distances for metal-poor 
stars are underestimated, while that for the metal-rich ones 
are slightly overestimated: this gives rise to a kinematic bias 
acting in opposite directions. 

5.5 Systematics of the approximation Adg = 0.7 

Fig-Slshows the discrepancy in heliocentric distance between 
the approximation Mg = 0.7 and our calibration on the two 
samples FS12 and Xll: the constant absolute magnitude ap- 
proximation tends to underestimate stars' distances by sev- 
eral kpc, the discrepancy is most severe for the most metal- 
poor objects. A direct assessment of the constant absolute 
magnitude approximation via the method in H2.ll measures 
a global distance underestimate of / = —0.13 ± 0.03 on Xll 
sample and of / = —0.11 ± 0.04 on FS12 sample, the error 
being worse for metal-poor objects. In a kinematic analysis 
for example, this systematics can distort the velocity distri- 
bution of the sample and the fact that the effect is stronger 
on metal-poor objects can introduce a false gradient in the 
kinematics of sub-samples of stars with different metallici- 
ties. 

We show that the method presented in H2.1l is sensitive 
enough to detect the metallicity and colour dependence of 
absolute magnitudes: for this purpose we bin both the spec- 
troscopic and the photometric samples separately in colour 
and metallicity to measure the expected dependences on 
each quantity while marginalising over the other. 

5.5.1 Colour dependence 

We divide Xll and FS12 into metal-poor, [Fe/H] G [-3, -2], 
and metal-rich stars, [Fe/H] G [—2,0]: each sub-sample is 
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Figure 4. Density distribution (p) of the difference in distance 
calibrations for BHB stars from SDSS DR9: As = s{Mg = 0.7) - 
s{Mg{{g - r)o, [Fe/H])). Solid lines refer to Xll sample, dotted 
ones to FS12. The assumption of a constant absolute magnitude 
underestimates stellar distances w.r.t. the colour and metallicity 
dependent estimate for Mg given by equation[5] The effect is more 
pronounced for the most metal poor objects (blue lines). 



then split in colour bins, the sizes of which are adjusted to 
provide enough statistics for each colour band, and therefore 
meaningful error bars 400 stars per bin). 

In Fig. [S]we plot / against [g — r)o: the red lines in the 
upper and lower panel are associated with the metal-rich 
stars in the Xll and FS12 samples respectively. The average 
relative distance error / oscillates with colour: the trend is 
most pronounced in the Xll sample. This behaviour reflects 
the dependence of their absolute magnitudes on colour. It is 
not matched by the metal-poor part of either of the samples: 
the blue lines in Fig. [S] show that both the photometrically 
selected metal-poor stars (lower panel) and the ones spec- 
troscopically selected (upper panel) exhibit an increase of 
/ with colour. The discrepancy in the colour-dependence 
at different metallicity is the footprint of the metallicity de- 
pendence of Mg. We remark that the metal-poor stars in the 
Xll sample (Fig. [S] upper panel, blue line) exhibit exactly 
the behaviour predicted in the previous section and are in- 
deed associated with the poorest fit by a constant absolute 
magnitude. Also, it is at these metallicities that our dis- 
tance calibration differs most from the one associated with 
Mg — 0.7: Fig. U clearly shows that our approximation pre- 
cisely corrects for the underestimate in distance caused by 
the assumption of a too faint constant magnitude. 

To assess the significance of these dependences w.r.t. the 
traditional level of uncertainty associated with the absolute 
magnitude of BHB stars (e.s.eA4^ = O.lS.lsirko et al.ll200"i: 
De Propris et al.ll2010l : iDeason et al.ll2011al : iNewberg et all 
20091 ) we need to relate the fractional distance error to the 



effective error in absolute magnitude, SMg = M'g—Mg. Given 
that 



Mg 



go -5 logj 



0.01 kpc 



(6) 



and s' — {1 + f)s, where s is the real heliocentric distance, 
then 
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Figure 5. Fractional average distance error as a function of colour 
at different metallicities under the assumption that Mg = 0.7 of 
the spectroscopic sample Xll (upper panel) and the photometric 
sample FS12 (lower panel). The green lines are the full sample, 
while the red lines are associated with metal-rich stars and the 
blue ones with the most metal-poor objects. The actual horizon- 
tal coordinate has been offset for each line for the error bars to 
be distinguishable. Every data point is associated with the centre 
of a colour bin: its size varies in order for every bin to contain 
a comparable amount of stars and therefore produce meaning- 
ful statistics. Grey dotted lines represent a n absolute magnitude 
uncertainty of ±0.18 ("e.g. ISirko et al.ll2004h . 



M'g 



go -5 logi 



0.01 kpc 



M,-51ogio(l + /) (7) 



so: £A/g ~ — 51ogjQ(l + /) . An absolute-magnitude 
uncertainty of ±0.18 will therefore correspond to / G 
[—0.08,0.09]. Fig. [5] shows that the fractional distance er- 
ror trends as a function of colour for different metallicities, 
significantly exceed this band (grey dotted lines). This im- 
plies that the effects of approximating the absolute magni- 
tude of BHB stars with Mg = 0.7 go beyond the generally 
acknowledged uncertainty of this approximation. 

Further, the effects of distance misestimates superpose 
as they ought to: in Fig. [Sj the green lines that represent the 
full samples are consistent with being averages of the blue 
and red lines associated with the metal-poor and metal-rich 
sub-samples respectively. Also, both samples confirm that 
the blue end experiences the worst fitting by a constant 
absolute- magnitude distance scale. The dramatic drop in the 
metal-rich side of the photometric sample (Fig. [5] bottom 
panel, red line) is however beyond predictable behavior and 



A new calibration for the Blue Horizontal Branch 7 



we therefore flag this end as possibly affected by biases out- 
side our control. Unfortunately, the low statistics associated 
with this sample do not allow a more detailed investigation 
of the problem. 

Fig. [6] shows the dependence of f on (g — r)o colour 
on both samples for our calibration and the constant ab- 
solute magnitude approximation. The performance of ((S)) 
on the spectroscopic sample Xll (blue solid line) is almost 
perfectly satisfactory, while FS12 (blue dotted line) appears 
slightly problematic at the blue end as expected. We cau- 
tiously avoid extensive comment on the photometric sample 
as regards colour dependence due to the low statistics avail- 
able for the blue end and the large uncertainties in the stellar 
parameters at that end. 



5.5.2 Metallicity dependence 

Fig. [7] shows the dependence of / on metallicity for both 
distance calibrations. |/| increases towards underestimates 
in distance (negative /) with decreasing metallicity in both 
samples (Fig. [T]). The discrepancy between the two distance 
calibrations is also affected by metallicity: still in Fig. [71 
we see that the offset between the blue lines (associated 
with the colour and metallicity dependent calibration) and 
the red lines (constant absolute magnitude) is largest at the 
metal-poor end for both the FS12 and Xll samples. 

The fact that the most metal-poor objects suffer the 
largest inaccuracy in distance (already visible in Fig. [5l up- 
per panel), and experience the largest discrepancy between 
the two distance calibrations explains why in Fig.[6]the offset 
between the two distance calibrations is significantly larger 
for the Xll sample w.r.t. the FS12 one: the former sample 
is more metal-poor than the latter and the different mixture 
of metallicities translates into a larger difference between 
Mg = 0.7 and our formula ©. 

Finally, we note that when we adopt the colour and 
metallicity-dependent distance calibration, / lies within the 
traditional uncertainty band associated with SMg = ±0.18 
(corresponding to f — —0.08 and / = 0.09 respectively) for 
both the spectroscopic and photometric sample and both 
w.r.t. colour and metallicity. In particular we tested the hy- 
pothesis that / is uncorrelated with metallicity when (O is 
adopted: a chi-squared test confirms that the discrepancy 
between the data and a constant trend is not enough to re- 
ject the null-hypothesis, 10-20 % significance level (s.l.) for 
Xll and 30-50 % s.l. for FS12. 

On the other hand the same hypothesis is clearly re- 
jected for the trends associated with Mg = 0.7 implying a 
statistically significant correlation between / and [Fe/H]: at 
the 1% s.l. for Xll and at the 5 - 10% s.l. for FS12. 



5.6 Improvement w.r.t. previous calibrations 

In Fig.[8]we show a comparison between the calibrations con- 
sidered above: the plot is an extension of Fig. [T] in the sense 
that it now shows the trend of Mg with {g — r)o colo ur at two 
differe n t metallicit i es als o for the calibrations of | Sirko et al.l 
(|2004l ): IXue et all l|201ll ) and lDeason et all l|2011bl ). We do 
not re-plot the isochrones to avoid crowding in the figure. 

We find that our calibration is the only one among the 
ones listed above not to suffer the "falling sky" pathology 
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Figure 6. Fractional average distance error as a function of colour 
for BHB stars from SDSS DR9. Blue lines refer to the colour and 
metallicity dependent distance calibration JSJl, red lines depict 
the constant Mg = 0.7. Solid lines and lighter colours are associ- 
ated with the spectroscopic sample Xll, dotted lines and darker 
colours with the photometric sample FS12. The discrepancy be- 
tween the two distance calibrations is larger for Xll due to the 
sample being more metal-poor than FS12 (see text). Grey dot- 
ted lines represent the quoted absolute magnitude uncertainty of 
±0.18. 
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Figure 7. Same as Fig.|6] but here we plot the fractional average 
distance error as a function of metallicity for the two samples. 



of SBA12 (see i |2.1[l . with the other calibrations showing 
the most severe effects on the metal-poor sub-sample, as 
expected from the previous assessment. In Fig. [9] we show 
the average correlation between W and sin I sin b cos b for 
the metal-poor ob jects in t he X l l sample for our calibra- 
tion, the one of Sir ko et all (|2004 ): IXu"e et al.l (|201ll ') and for 
the constant magn itude approximation. The performance of 
the calibratio n oflDeason et al, (1201 lb[) is almo s t iden tical 
to the ones of lSirko et al.l l|2004) and lXue etall (120111 ) and 
we avoid plotting it in order not to crowd the figure. We 
note that the offset from zero in the mean W disappears 
accounting for a vlos aberration of ~ 10 — 15kms^^. While 
this is not a proof that SEGUE viob suffer from a system- 
atic offset, the coincidence requires further investigation: we 
reserve this for future work. 
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Figure 8. Same as for Fig. [T] but including the calibrations of 
Sirko et al. (2004), Xue et al. (2011) and Deason et al. (2011b). 
The calibration of Deason et al. (2011b) is independent of metal- 
licity and hence shown only once. 
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Figure 9. Mean correlation between the W component of velocity 
and the angular term sin i sin b cos b for metal-poor stars ( [Fe /H] < 
—2) in the Xll sample. Solid lines are the mean trends associated 
with the different calibrations: green for our eq. (5] red for the 
calibration of Sirko et al. (2004), Xue et al. (2011) and blue for 
the constant absolute magnitude approximation. With the same 
colour code, dashed lines represent the Icr contours. 

6 CONCLUSIONS 

We have proposed and validated a very simple analytic 
formula that properly describes the colour and metallic- 
ity dependence of the absolute magnitudes of BHB stars. 
Our calibration both provides a user-friendly tool to com- 
pute the absolute magnitude of BHB stars from observ- 
ables in the ugriz-system and achieves the best performance 
on field BHB halo stars w.r.t. other available calibrations. 
Given the agreement between our formula and the BASTI 
isochrones, we conclude that these isochrones successfully 
predict ugriz colours for the blue part of the HB. The Dart- 
mouth isochrones are fainter than the BASTI models by 
~ 0.1 — 0.4 mag and therefore lead to underestimate the dis- 
tances of field BHB stars. 

The quality of proper motions in the SDSS is sufficient 
to assess the average distance scale. Even though the large 



distances covered by the two samples imply large random er- 
rors in kinematics from proper motions, our considerations 
are based on averages across the sky. The linear estimator we 
adopted from SBA12 is therefore not affected by the inac- 
curate knowledge of the size of proper- motion errors, while 
it is mildly vulnerable to p roper moti o n sys tematics. We 
corrected for them following ISchonrichI l|2012l ) and showed 
that this correction leads to a stretch of the distance scale 
of about 5%. 

Another issue is the unknown level of contamination 
by blue stragglers and main-sequence stars. This contam- 
ination will bias our distance estimates towards signalling 
an overestimate. Hence our formula providing formally un- 
biased distance estimates may still be rather on the faint 
side despite being brighter than the alternative calibrations. 
Our magnitude approximation may be too flat at the blue 
end of the horizontal branch, where it declines towards the 
blue hook. We have actively chosen this mild deviation in 
the light of reddening and photometric uncertainties. 

The fact that the absolute magnitude of BHB stars de- 
pends on both colour and metallicity implies that any cali- 
bration that does not account for this dependence is a pri- 
ori biased. We have shown that impact of a constant mag- 
nitude approximation on distances exceeds the 10% level 
(/ G [-0.38,0.04] or equivalently AMg G [-0.09,1.04], see 
Fig. [5l [7]) and that a calibration which is dependent on 
colour, but not on metallicity (e.g. iDeason et al.ll2011bh m- 
troduces a systematics in distance (of order of 0.23 mag or 
respectively 10%) which correlates with metallicity. This ef- 
fect is of particular concern in kinematics studies for it can 
produce a false kinematics gradient across populations with 
different metallicity. 

The availab le calibrations that depend on both co lour 
and metallicity l|Sirko et al.ll2004l : IXue et al.ll2008l . l201ll ) re- 
quire interpolation of theoretical isochrones for each magni- 
tude estimate and they were validated on a handful clusters, 
which themselves have a distance uncertainty. We showed 
that on field BHB stars they significantly under-perform our 
suggested calibration (on metal-poor stars / can be as high 
as 0.17) and actually fail the "falling sky" test of SBA12. 

The simple calibration formula we presented in equa- 
tion [5] of this paper proves to be more accurate than any 
of the above on field BHB stars, achieving an accuracy that 
matches the threshold of the noise determined by the sample 
size, being the only one that passes the "falling sky" test of 
SBA12. Larger information content of the data (e.g. better 
metallicity estimates or a new dimension like Helium con- 
tent) and better statistics would be required to improve the 
accuracy of our formula and assess its performance in the 
blue-hook region. 

In general, the method we use to assess distance cali- 
brations can be applied to any filter set and sample of field 
stars in the sky. The only requirements are that the sam- 
ple comprises a sufficiently large number of objects with 
acceptable proper motion quality, is sufficiently extended in 
Galactic coordinates and especially is not associated with a 
strong kinematic selection or is dominated by a low number 
of streams. 
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